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Abstract 

Alterations in gut microbiome composition have an emerging role in health and disease including brain function and 
behavior. Short chain fatty acids (SCFA) like propionic (PPA), and butyric acid (BA), which are present in diet and are 
fermentation products of many gastrointestinal bacteria, are showing increasing importance in host health, but also may be 
environmental contributors in neurodevelopmental disorders including autism spectrum disorders (ASD). Further to this we 
have shown SCFA administration to rodents over a variety of routes (intracerebroventricular, subcutaneous, intraperitoneal) 
or developmental time periods can elicit behavioral, electrophysiological, neuropathological and biochemical effects 
consistent with findings in ASD patients. SCFA are capable of altering host gene expression, partly due to their histone 
deacetylase inhibitor activity. We have previously shown BA can regulate tyrosine hydroxylase (TH) mRNA levels in a PCI 2 
cell model. Since monoamine concentration is known to be elevated in the brain and blood of ASD patients and in many 
ASD animal models, we hypothesized that SCFA may directly influence brain monoaminergic pathways. When PCI 2 cells 
were transiently transfected with plasmids having a luciferase reporter gene under the control of the TH promoter, PPA was 
found to induce reporter gene activity over a wide concentration range. CREB transcription factor{s) was necessary for the 
transcriptional activation of TH gene by PPA. At lower concentrations PPA also caused accumulation of TH mRNA and 
protein, indicative of increased cell capacity to produce catecholamines. PPA and BA induced broad alterations in gene 
expression including neurotransmitter systems, neuronal cell adhesion molecules, inflammation, oxidative stress, lipid 
metabolism and mitochondrial function, all of which have been implicated in ASD. In conclusion, our data are consistent 
with a molecular mechanism through which gut related environmental signals such as increased levels of SCFA's can 
epigenetically modulate cell function further supporting their role as environmental contributors to ASD. 
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Introduction 

The gut microbiota - the diverse range of symbiotic gut bacteria 
and other microorganisms is involved in the regulation of multiple 
host metabolic pathways in both health and disease [1,2]. There is 
increasing evidence this microbial ecosystem, which outnumber 
host cells 100 to one, act as a functional "organ", playing a major 
regulatory role in gut-brain communication, brain function and 
even behavior [3,4,5,6,7]. The mutually beneficial relationship 
between the host and gut microorganisms arises in part from 
SCFAs which are produced from bacterial fermentation of some 
proteins and dietary fiber, the most abundant of which are acetate, 
BA and PPA [8] . These SCFA serve local functions in phenotypic 
reprogramming of colonic epithelial cells, as the principal energy 



substrate for epithelial cells, as tumor suppressor agents, in 
apoptotic cell death, and in gene regulation of anti-inflammatory 
processes both in vitro and in vivo [9,10,11,12,13,14,15]. 

The majority of SCFA are absorbed from the gut, transported 
via the portal vein and metabolized in the liver before reaching the 
systemic circulation. Thus, hepatic clearance was thought to 
impede achieving levels sufficiently high to affect systemic 
functions or metabolic and regulatory pathways. However, the 
distal colon bypasses the portal circulation enabling systemic 
access. Despite the difficulty in their measurement due to 
intracellular concentration and rapid metabolism, there is growing 
evidence that the potential systemic effects of SCFA (especially 
PPA and BA) on physiology and pathology may have been 
underestimated (rev. in [8,3]). 
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SCFA have a number of direct efiects on gastrointestinal 
physiology. Along with acetate and butyrate, PPA is known to 
reduce gastric motility and increase the frequency of contractions, 
presumably via a reflex that involves direct contact of these short 
chain fatty acids with the terminal ileum [16]. In addition, PPA 
increases contraction of colonic smooth muse Ic [17], dilates 
colonic arteries [18], activates mast cells [19] and increases the 
release of serotonin from gut enterochromaffin cells [20,21]. In 
spite of the multiple beneficial effects of SCFA on host 
gastrointestinal activity, excessive quantities of PPA have been 
reported in gingival inflammation [22], acne [23], irritable bowel 
syndrome [24] and in the neurometabolic disorder propionic 
acidemia [25]. In this heterogeneous inborn error of metabolism, 
which may be underreported [26], accumulation of PPA is 
associated with developmental delay, seizure and extrapyramidal 
findings, often accompanied by gastrointestinal symptoms [25,27], 
confirming that a homeostatic balance of SCFA may be necessary. 
Furthermore PPA and related SCFA have broad effects on 
nervous system physiology, including activation of specific G 
protein coupled receptors (GPCR), neurotransmitter synthesis and 
release. intraix'Uular pll/calcium gating, mitochondrial function, 
lipid metabolism, immune function, gap junction gating and gene 
expression (see [3]). 

Recent evidence suggests potential, but unproven, links between 
dietary, metabolic, immune, infective, and gastrointestinal factors 
and ASDs. Although inheritable factors, mostly implicated in 
synaptic transmission, have been traditionally studied in ASD [28], 
the fact that 1) known genetic factors thus far account for only 10- 
20% of cases, 2) there is less than 100% concordance in identical 
twins, and 3) there is a growing prevalence in the condition, 
collectively suggest an important role for environmental factors 
which act on the underlying genetic sensitivities [29] . In particular 
hospitalization, early infections and associated antibiotic exposure 
[30], which are risk factors for ASD, may alter the developing gut 
microbiome [2] . Increased mean levels of PPA in stool of ASD 
children have been shown [31] although contrasting findings were 
also reported in ASD participants with higher usage of probiotics 
and fish oil consumption [32]. Given that PPA is a key 
fermentation product of ASD-associated antibiotic resistant 
bacteria {Clostridia, Bacteriodetes, Desulfovibrio) [33] and modu- 
lates many ASD related biochemical processes, we have proposed 
that SCFAs represent a group of host microbiome metabolites that 
are plausibly linked to ASDs and can induce widespread effects on 
gut, brain, and behavior [3,4]. 

Furtiier to this, we [34,35,36,37,38,39,40] and otiiers 
[41,27,42,43,44] have shown that short-term central nervous 
system (intracerebroventricular) and peripheral administration 
(intraperitoneal, subi:utaneous or oral gavage) of PPA, and to a 
lesser extent, other SCFAs at various developmental time periods 
in rodents induce ASD-like abnormal motor movements, repet- 
itive interests, electrographic changes, cognitive deficits, persever- 
ation, and impaired social interactions. The brain tissue of PPA- 
treated rats also shows a number of ASD-linked neurochemical 
changes, including neurotransmitter alterations, innate neuroin- 
flammation, increased oxidative stress, redox changes, glutathione 
depletion, and altered phosphoUpid/acylcarnitine profiles. This 
novel model also has shown predictive value for potential 
metabolic biomarkers in a large cohort of ASD patients [45]. 

One potential key mechanism where the metabolic products of 
an altered gut microbiome may contribute to ASD pathophysi- 
ology is via the alteration of gene expression associated with ASD 
mutations or ASD implicated genetic pathways [3,46,7]. Interest- 
ingly, PPA, related SCFA and their derivatives are known 



modulators of gene expression via their histone deacetylase 
inhibitor (HDACI) activity [47,48,49,50,51,52,53]. 

The rat pheochromocytoma (PCI 2) cell line is an extensively 
used in vitro cell system to examine molecular biological processes 
in neurobiology [54] . We have used this PC 1 2 line to examine the 
effects of SCFA, principally BA, and their derivatives on gene 
expression [55,56,57,58]. Of note, our results underscored at least 
3 major mechanisms by which BA can regulate TH gene 
expression: i) modulation of gene transcription by chromatin 
remodeling, ii) by activation of transcription factors via different 
signahng cascades (including Ca2+/cAMP mediated activation of 
CREB [59]) and involving induction of transcription via an 
upstream 5' regulatory element (BRE; GCCTGG at -509 to - 
504 of the rat TH promoter [60] or iii) by affecting the stability of 
TH mRNA (e.g. via a butyrate response factor (BRF) acting at a 3' 
untranslated AU-rich regions of mammalian mRNAs 
[61,62,63,64]). Many of the affected genomic pathways are 
involved in catecholamine synthesis, which have been implicated 
in ASD [65,66]. Moreover, CREB, a key factor in neurodevel- 
opment, learning and memory [67], is a key determinant of 
catecholamine synthesis in PC 12 cells, and shows increased CREB 
immunoreactivity in brains of PPA treated rats (our animal model 
of ASD, [35]). Furthermore, the anti-seizure/mood stabilizing 
drug valproic acid, a known prenatal risk factor for ASD, which 
produces an acceptable animcd model for the condition, is 
structurally and pharmacologically similar to PPA, including 
HDACI properties [68,69,70] and produces similar effects as BA 
in PC 12 cells [58]. 

Since most research on the effects of SCFA on gene expression 
is limited to BA and not PPA, in the present study we used rat 
PC 1 2 cells as an in vitro system to extend our observations on the 
epigenetic effects of SCFA. Microarray technology was used to 
compare global changes in gene expression profiles following 
exposure to the structurally related SCFA PPA and BA. 
Furthermore, we sought to determine if the expression of these 
PPA dependent genes was related to canonical biological pathways 
implicated in ASD. 

Methods 

Cells/Transfection 

PC 12 cells ([71], rat pheochromocytoma of sympathoadrenal 
origin) were used as a model to deUneate the molecular effects of 
propionate. They were cultured in DMEM media supplemented 
with 1 0 % horse serum, 5 % fetal bovine serum and antibiotics in a 
humidified atmosphere and lO'Vo CO2 as described earlier [64]. 
Sodium propionate (Sigma, St. Louis, MO) was added at the 
indicated concentrations at nearly 50% confluence. To evaluate 
the effect of PPA on TH gene transcription we used a transient 
transfection approach. Briefly, plasmid constructs with a rat TH 
promoter (-7737+27 bp, [72] driving the expression of firefly 
luciferase reporter gene were electroporated into PC 12 cells at 
300 V and 500 [iF with a total of 50 ng/plate of plasmid DNA 
(lO^-lO' ceUs). CeUs treated with SCFA continued to receive 
media supplemented with SCFA after electroporation. After 
incubation for an additional 24 hours, ceUs were harvested and 
crude cell lysates were prepared. Protein concentrations in total 
ceU lysates were determined [73] and luciferase reporter gene 
activities were measured by luminometry as recommended by 
Promega (Madison, WI). AU averages were made from at least six 
independent assays. The results were given as relative light units 
per second per microgram of total protein (RLU/ s/ |xg). The 
luciferase activities were expressed as fold stimulation relative to 
ceUs co-transfected with the same constructs but not treated with 
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SCFA. The values shown are means ± SEM from three 
experiments. 

Wild type or dominant negative CREB expression vectors were 
purchased from BD Biosciences (Palo Alto, CA). AU plasmid DNA 
constructs were purified by chromatography on Qiagen columns 
(Santa Clarita, CA). 

Isolation of RNA and Northern blot analysis 

To compare the effects of structurally related SCFA on 
endogenous TH gene expression, PC 12 cells were treated with 
low concentrations of SCFA (PPA, BA, and VP A) or with vehicle 
for 48 hours. Total RNA was isolated from individual petri dishes 
using RNazol according to the manufacturer's protocol (Tel-Test, 
Inc., Friends-Wood, TX). Northern blot analysis was performed as 
described previously [64]. After transfer to Gene Screen Plus 
membranes, the filters were hybridized consecutively to labeled rat 
cDNA probe for TH and a probe for 1 8S ribosomal RNA. Filters 
were washed under proper stringency and exposed for autoradi- 
ography, using Kodak Bio Max film (Rochester, NY). The 
autoradiographs were scanned and quantified by Bio-Rad 
Quantity One software and the abundance of TH mRNA was 
expressed relative to 18S ribosomal RNA levels. The results were 
presented as fold change compared with the control (vehicle 
treated) group on the same Northern blot. 

Microarrays 

PC 12 ceUs were treated with 1 mM PPA, 10 mM PPA, 1 mM 
BA, 6 mM BA or vehicle for 48 hrs. Total RNA samples from 
each experimental group were pooled together and after DNase 
treatment purified using RNeasy Mini Kit (Qiagen Sciences, MD) 
as recommended by Afifymetrrx. The cDNA synthesis and 
microarray analyses were performed at Keck Affymetrix Gene- 
Chip Resource at Yale, New Haven, CT (NIH Neuroscience 
Microarray Consortium). The chip data was analyzed using 
Affymetrix GeneChip Operating Software (GCOS) version 1.2.1. 
Genes showing altered expression with fold change values of ^2.0 
were exported for functional annotation to MetaCore by 
GeneGO, Inc. The expression levels of selected genes (TH, 
DBH, NPY, ChrA, PENK, and GTPCH) were validated through 
Northern blot analyses of the same RNA samples (Table 1). Raw 
and quantile-normalized microarray data and an associated 
project metadata file are available through the NCBI-GEO 
repository (GSE56516). 



Western blot analysis 

Total protein extracts were prepared from each petri dish of 
control (vehicle treated) and SCFA- treated PC 12 cells for the 
indicated periods of time (see the figure legend). Proteins were 
separated on 10% SDS-PAGE, electroblotted onto a nitrocellulose 
membrane (BioRad; Hercules, CA) and incubated with TH 
antibody (1: 4000, Imgenex, San Diego, CA) overnight. After 
incubation with secondary antibody (Goat Anti-Rabbit IgG, from 
Pierce, Rockford IL; diluted 1: 40000) the immune reaction was 
visualized by enhanced chemiluminescent substrate from Pierce, 
utilizing a horseradish peroxidase label and Kodak XAR-5 film, as 
described by the manufacturer. To confirm equal loading, blots 
were re-probed with primary antibody for the house keeping 
protein fi-actin (Monoclonal Anti-fi-Actin antibody, Sigma, St. 
Louis, Mo). The blots were exposed to autoradiography and the x- 
ray films were scanned and quantified with Bio-Rad Quantity One 
software. The ratio TH/fi-actin immunoreactivity were calculated 
for each sample and the results were presented as fold induction 
compared to the corresponding control group on the same 
Western blot. 

Statistical analysis 

Data were expressed as mean ± SEM and normalized to the 
values in the control, taken as 1.0. Differences between the 
experimental groups from at least three independent experiments 
were evaluated by performing ANOVA followed by the Fisher's 
least significant difference test for experiments with more than two 
groups. A level of pSO.05 was accepted as statistically significant. 

Results 

1. PPA induces TH gene transcription in a CREB- 
dependant fashion 

We have shown before [63,64,60,58] that HDACI like BA and 
VPA can also activate transcription of the TH gene via upstream 
regulatory elements in PC 12 cells. To test whether PPA has similar 
effects we performed transient transfection experiments following 
the same protocol. Briefly, plasmids carrying the luciferase 
reporter gene under the control of the rat TH promoter (— 
773/-H27 bp) were electroporated into PC 12 cells pre-treated for 
one day with vehicle or with SCFA. Reporter gene activity was 
measured in control and SCFA-treated culture samples 24 hours 
after transfection. PPA increased luciferase activity in a concen- 
tration-dependent manner consistent with activation of TH 
promoter (Figure lA). Furthermore, introducing a single point 
mutation into the canonical CRE motif (a G-to-A alteration at 



Table 1. Effect of PPA and BA on neurotransmitter related gene expression In PCI 2 cells. 





Gene/fold change 


1 mM BA 


6 mM BA 


1 mM PPA 


10 mM PPA 


TH 


t 2.8±0.21 


i 6.4±0.12 


t 2.17±0.21 


i2.2±0.35 


DBH 


t 1.4±0.17 


I 14.5±0.89 


t 1.2 ±0.23 


J, 1.8±0.13 


NPY 


t 2.78±0.35 


t 5.2±0.64 


1.02 ±0.28 


t 1.3±0.35 


ChrA 


t 1.8±0.15 


I 3.0 ±0.48 


t 1.15±0.11 


J, 0.85±0.12 


PENK 


t 4.65 ±0.39 


t 2.46±0.22 


t 2.7±0.25 


t 1.95±0.25 


GTPCH 


1.02±0.09 


0.98±0.15 


t 1.45 ±0.23 


t 2.67±0.27 



PCI 2 cells were treated witfi 1 mM PPA, 10 mM PPA, 1 mM BA, 6 mM BA or vehicle for 48 hrs. Total RNA samples were isolated from individual petri dishes (n>6 per 
experimental group). Nortfiern blot analysis was performed as described in Metfiods section, using probes specific for rat genes, encoding tyrosine hydroxylase (TH), 
dopamine beta-hydroxylase (DBH), neuropeptide Y (NPY), Chromogranin A (ChrA), pre-proenkephalin (PENK) and GTP cyclohydrolase {GTPCH}. The results are 
presented as mean ± SEM value relative to the mRNA levels in control, vehicle treated cells. 183 rRNA was used as a house keeping gene for normalization of data. 
doi:l 0.1 371/journal.pone.Ol 03740.t001 
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Figure 1 . PPA activates the transcription of TH gene in CREB-dependant manner. Plasmid constructs with wild type (WT) rat TH promoter 
(—7737+27 bp) driving the expression of luciferase reporter gene or having a G-to A transition mutation in the CRE promoter element (mCRE) were 
electroporated in PCI 2 cells pre-treated for one day with increasing concentrations of PPA (0.1 to 10 mM, Fig. 1A)). The same doses of PPA were 
given immediately after electroporation for additional 24 hours, when changes in reporter gene activity were determined. In the second set of 
experiments (Figure IB), PCI 2 cells were co-transfected with combination of wt TH promoter construct and plasmids expressing wild type or 
dominant negative CREB mutants (binding-defective K-CREB or phosphorylation-defective 133CREB). To half of the cultures, 10 mM PPA was given 1 
day before and immediately after electroporation for additional 24 h (black bars). Controls received vehicle 24 hrs before and immediately after 
electroporation (open bars). The values shown are means ± SEM from three experiments. The luciferase activities (RLU/s/|j,g protein) are expressed as 
fold stimulation relative to control cells co-transfected with the same plasmids but not treated with PPA. *P<0.05 (PPA vs. vehicle treated; **P<0.002 
(wt vs. mutant constructs). 
doi:1 0.1 371/journai.pone.01 03740.g001 



position —41, mCRE) of the TH promoter caused more than 85% 
inhibition of the PPA-induced reporter gene activity. We had 
similar results obtained previously for BA and VPA [60,58], 
suggesting involvement of Ca^^^cAMP-mediated signaling path- 
ways in the transcriptional effects of structurally related SCFA. 

To test this possibility, we also evaluated the role of the general 
transcription factor CREB in PPA-induced activation of TH gene 
promoter. In these experiments, PCI 2 cells were transiently 
transfected with plasmids carrying the luciferase reporter gene 
under the control of rat TH promoter in combination with either 
wild type CREB (positive control); dominant negative CREB 
expression vector(s) or unrelated DNA (pCMV-CAT expression 
vector) as we described before [56]. PPA was given one day before 
and immediately after the transfection. Over-expression of non- 
functional CREB mutants unable to bind to their cognate DNA 
enhancer (K-CREB) or to get transactivated in response to 
activation of different signaling cascades (133CREB) significandy 
reduced the PPA-triggered increases in reporter gene activity 
(Fig. IB). These results indicate that CREB or a related protein(s) 
and signaling pathways converging on them are involved in 
mediating the transcriptional effects of PPA on the TH gene 
promoter. 



2. PPA can alter steady state TH mRNA and TH protein 
levels 

Our previous data revealed that exposure of PC 1 2 cells to low 
concentrations of butyrate resulted in accumulation of endogenous 
TH mRNA (see [64,63]). To compare the effect of structurally- 
related SCFA on TH mRNA levels, PC 12 cells were treated 
1 mM BA, 1 mM PPA or 0.1 mM VPA for 48 hours. All 
treatments resulted in a statistically significant increase in TH 
mRNA levels (Fig. 2A). 

To determine whether SCFA-induced alterations in steady state 
TH mRNA levels were reflected in changes in TH protein, we 
treated PC 1 2 ceUs with 1 mM PPA or BA for 24, 48 and 72 hours. 
Controls received only vehicle. Exposure to SCFA resulted in 
gradual increases in the relative immunoreactive TH protein 
levels, suggesting that the capacity of the cells to produce 
catecholamines is elevated (Figure 2B). 

3. Structurally related SCFA regulate similar genetic 
networks in PCI 2 cell model 

To identify other potential cellular targets for SCFAs, we used 
microarray technology allowing genome-wide simultaneous mea- 
suring of changes in gene expression. PC 1 2 cells were exposed to 
SCFA (1 mM BA, 6 mM BA, 1 mM PPA and 10 mM PPA), or 
vehicle for 48 hours. Total RNA was isolated and subjected to 
microarray analysis using Affymetrix GeneChip Rat Genome 230 



PLOS ONE I www.plosone.org 



4 



August 2014 | Volume 9 | Issue 8 | e103740 



Epigenetic Regulation of Neuronal Genes by Short Chain Fatty Acids 




C PPA BA VPA 0 24 48 72 hrs 

Figure 2. PPA induces accumuiation of TH mRNA and Tl-i protein in vitro. PC12 cells were treated with SCFA: propionic acid (PPA, 1 mM), 
butyric acid (BA, 1 mM) or valproic acid (VPA, 0.1 mM). At indicated times (48 hrs) total RNA or total cell lysates (0, 24, 48 and 72 hrs) were prepared 
and subjected to Northern (A) or Western blot analyses (B, PPA and BA). The data on Fig. 2A are presented as mean ± SEM value relative to the mRNA 
levels in control, vehicle treated cells. 18S rRNA was used as a house keeping gene for normalization of data. *P<0.05 SCFA vs. vehicle treated. For 
Western blotting (23) cells were lysed after indicated incubation times and equal amounts of proteins were subjected to analyses using TH-specific 
antibody as described In Experimental procedures. Membranes were re-probed with p-actin specific antibody as loading control. The 
autoradiographs were analyzed by densitometry. Values are presented as the ratio of TH to that of beta-actin as fold difference from control, 
vehicle- treated cells. N = 6-8 independent samples per group, and the experiment was repeated twice with similar results [63]. 
doi:l 0.1 371/journal.pone.0103740.g002 



2.0 microarray as described in Methods. Data was ranked 
according to significant up-regulation or down-regulation of 
twofold or more in SCFA-treated groups compared to their 
respective control groups (*P<0.05). 

Only few affected genes (79) fit the criteria when cells were 
treated with low dose of 1 mM PPA. The expression of a 
significantly higher number of genes was altered following 
exposure to 10 mM PPA (1587). At both doses, BA was more 
eflFective than PPA as it was evident by the number of dififerentially 
expressed transcripts and correlated with the higher potency of B A 
to induce histone hyperacetylation and alter global gene expres- 
sion. For example, a total of 42 1 genes were significantly altered 
following treatment with 1 mM BA, and 3148 for 6 mM BA 
group. We report here the results from MetaCore enrichment 
analysis performed for 10 mM PA and 6 mM BA groups. 

Comparison analysis revealed 1599 differentially expressed 
genes unique for the BA experimental group and 355 for the 
PPA group. A large number of genes (1010) were identified as 
common for both, BA and PPA groups (Fig. 3, complete lists 
provided in Table SI, S2, S3). Examples of the highest magnitude 
affected common genes are given on Table 2. In most cases the 
extent of induction/ reduction of gene expression triggered by BA 
was more pronounced compared to PPA with the exception of 
only 16 genes (Table 3). There was also a small subset of genes (21) 
that change in opposite direction for PPA and BA (Table 4). 

4. Specifically Affected Gene Systems 

Differentially expressed common genes included but were not 
limited to the following: those involved in the biosynthesis and 
degradation of dopamine, norepinephrine and serotonin (i.e. TH, 
DBH, DDC, TPH, GTPCH, COMT, MAOA), mitochondrial/ 
fatty acid metabolism (i.e. 3-ketoacyl-CoA thiolase, ACAA2; acetyl 
CoA acetyltransferaseACAT2; long chain fatty acid CoA ligase 5- 
ACSL5; 3-hydroxybutyrate dehydrogenase type 2, BDH2; Glyc- 



erol-3-phosphate dehydrogenase, GPDl; hormone sensitive lipase, 
LIPE), oxidative stress (GPX3-glutathione peroxidase; APOE- 
apolipoprotein E; MT3- metallothionein-3; DUSPl dual specific- 
ity protein phosphatase 1), gap junc:tion proteins, signaling 
pathways (cAMP/Ca2+, MAP kinase) receptors (nicotinic, 
GABA), see Tables 5-8. The complete Kst of commonly affected 
by SCFA genes is available in Table S3. For a few selected genes, 
the expression patterns were validated by Northern blot analyses 
(Table 1). 

5. Gene Networks 

The top five most relevant gene networks, identified using the 
Analyze Networks (AN) algorithm with default settings (MetaCore, 
[74]) are listed in Table 9. Processes involved in the cell cycle and 
regulation of the progression through the cell cycle are the most 
affected, consistent with the ability of HDACI to induce growth 
arrest, cell differentiation, and chromatin remodeling in a variety 
of ceU fines [47]. 

In order to explore the biological significance of SCFA- 
mediated gene expression changes, we analyzed their Gene 
Ontology (GO) classification using MetaCore (Fig. 4). The most 
highly represented processes (sorting method "similarity by") 
include: synaptic transmission, regulation of G-protein coupled 
receptor signaling, response to stimulus, transmission of nerve 
impulse, cardiac muscle development, regulation of NE secretion, 
NE secretion. 

Discussion 

Emerging evidence hnks SCFA from diet, cellular metabolism 
or pharmacological treatments with regulation of neuronal gene 

expression, attendant neuronal function and even behavior 
[34,35,36,37,38,45,39,41,42,27,43,44,64,58,63,75], rev. in [3]. In 
this study, we demonstrate that PPA can induce TH gene 
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# genes 



1599 




unique common 

355 



1010 



*p< 0.001 



Unique for PPA 
Unique for BA 



Figure 3. Differential gene expression profiles evoiced by SCFA in PC12 ceils. Total RNA was isolated from PC12 cells following 48 hrs 
exposures to BA (6 mM), PPA (10 mM) or vehicle and subjected to microarray analysis. Genes that changed two-fold or greater in SCFA-treated 
groups (*p<0.05) compared to the control (vehicle treated PCI 2 cells) group were imported into MetaCore GeneGo for comparison analyses. 
Differentially expressed unique (1599 for BA and 355 for PPA experimental group) and common for both groups genes (1010) were identified. Fisher's 
exact test p-value was <2.2e-16. 
doi:1 0.1 371/journal.pone.01 03740.g003 



transcription over a wide dose range (0.1 niM to 10 mM) in PC 12 
cells, similar to the elTects of BA [63,75]. The general transcription 
factor CREB (or CREB-related proteins) and signaling pathways 
converging on it are involved in the transcriptional regulation of 
the TH gene by PPA or BA. At lower concentrations, PPA - 
triggered changes in TH gene transcription are reflected in 
accumulation of TH mRNA and TH protein, which are markers 
of increased catecholamine production. Genome-wide expression 
profiling following exposure to SCFA enabled us to identify large 
number of affected genes and networks, common for both PPA 



and BA, including catecholaminergic pathways and GO processes 
implicated in ASD. 

The short chain fatty acid BA is a well known HDACI shown to 
induce differentiation, growth arrest and chromatin remodeling in 
a variety of cell lines, associated with specific-gene transcription 
regulation properties [47,48,49,50]. The effects of BA are similar 
to its structurally related branched short chain fatty acid with a 
known causal role in ASD VPA, [51,52] which produces an 
acceptable animal model for the condition [68,70,69,76,77]. PPA 
also exhibit HDACI properties in non-neuronal cell lines, 



Table 2. Comparison analysis of differentially expressed genes triggered by PPA and BA in PCI 2 cells: Top 10 most affected 
common genes. 



Symbol 


Gene 


PPA 


BA 


COL5A3 


Collagen AlphaB (V) chain 


t 644.5 


t 1 048.6 


COX6B2 


Cytochrome c oxidase subunit VIb 


t 6.634 


t 156 


CA5B 


Carbonic anhydrase 5B, mitochondrial 


t 8.381 


t 91.615 


CD40 


Tumor necrosis factor receptor superfamily - member 5 


t 33.828 


t 90.103 


ICAM1 


Intercellular adhesion molecule 1 


t 3.287 


t 89.463 


HPD 


4-hydroxyphenylpyruvate dioxygenase 


t 2.821 


t 82.462 


PPP1R1B 


Protein phosphatase 1 regulatory subunit IB 


t 2.158 


\ 54.368 


SH3KBP1 


SH3 domain-containing kinase binding protein 1 


t 2.8 


t 40.719 


SELENBP1 


Selenium-binding proteinl 


t 14577 


t 25.455 


SPHK1 


Sphingosine kinase 1 


t 5.177 


t 19.341 


TTK 


Dual specificity protein kinase TTK 


i 92.8 


i 189.1 


TOP2A 


DNA topoisomerase 2 alpha 


i 16.8 


i 178.7 


NDC80 


Kinetochore protein HEC1 


i 541 


i 170.2 


RRM2 


Ribonucleoslde -diphosphate reductase 


i 12.9 


i 160.9 


ALDOA 


Fructose-bisphosphate aldolase A 


i 22.88 


i 142.7 


CCDC95 


Coiled coil domain containing protein 95 


i 22.8 


i 142.7 


KIAAOIOI 


PCNA-associated factor 


i 22.88 


i 142.7 


GTSE1 


G2 and S-phase expressed protein 1 


i 7.6 


i 122.9 


MKI67 


Antigen KI-67 


i 16.631 


i 85.181 


ECT2 


Protein ECT2 


i 15.863 


i 82 



The symbols ] and [ indicate up- and down-regulated genes, respectively. Numbers following the f or j symbols represent the fold-change for the gene 
expression level. Sorting was done according to the magnitude of the signal (fold change) in BA group. The top 10 for each, upregulated and down regulated genes are 
shown. For more details, see Table S3. 
doi:1 0.1 371/journal.pone.Ol 03740.t002 
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Table 3. Examples of common genes which were affected more by PPA. 



Symbol Gene PPA BA 

CD7 T-cell antigen CD7 1 33 ] ^8.3 

FBP2 Fructose-l,6-biphosphatase isozyme 2 1 27.96 f^-lS 

STYXLI Serine/threonine/tyrosine-interacting like protein 1 ^12.9 f'9.2 

DYSF Dysferlin t 12.7 f 4.1 

CBS Cystatiiionine beta-synthase f 10.6 12.25 

TPPP3 Tubulin polymerization-promoting protein 3 f 6.82 1 4.2 

DRD1IP Dl dopamine receptor-interacting protein calcyon ^6.1 fS-l 

SEMASA Semaphorin 6A 1 6.1 1 2.7 

CGN Cingulin 1 5.98 12.26 

TUFT1 Tuftelin 1 5.98 12.2 

FOS Proto oncogene protein c-fos 1 5.86 f 2.2 

TCEA3 Transcription elongation factor A protein 3 1 5.82 f 3.77 

RKAG2 5'-AMP-activated protein kinase gamma 2 i 13.6 J, 2.1 

BARD1 BRCAl -associated RING domain protein i 7.8 i 4.4 

CCNE2 Gl/S specific cyclin-E2 i 4.48 i 2.8 

SYCP3 Synaptonemal complex protein 3 i 3.98 i 2.2 



Sorting was done according the magnitude of the signal (most affected, fold change) in PPA group for both upregulated { f ) and down regulated genes ( I ). 
doi:l 0.1 371 /journal.pone.Ol 03740.t003 

although it is less potent than BA [53]. Here using a PC 12 cell phenylbutyrate [75,58]. This induction was dependent on the 
model we show that PPA activates TH promoter-driven reporter binding ability and transactivation status of the general transcrip- 
gene expression in a concentration-dependent manner (Fig. lA), tion factor CREB - a stimulus induced transcription factor, 
similar to other HDACI like BA, VPA, trichostatin A and phosphorylated at Serl33 by diverse array of protein kinases 

Table 4. Common genes oppositionally affected by PPA and BA. 



Symbol Gene PPA BA 

SPCS3 Signal peptidase complex subunit 3 i 2.9 f 307.5 

GASe Growth arrest specific protein 6 1 6.6 f 20.48 

RGS2 Regulator of G-protein signaling 2 i 2.4 110.4 

HLA-G HLA class 1 histocompatibility antigen alpha G i 2.4 1 6.71 

NCAM1 Neural cell adhesion molecule 1 i 3.0 f 4.4 

HABP4 Intracellular hyaluronan binding protein 4 i 2.2 1 4.03 

BRD2 Bromodomain- containing protein 2 ^ 2.5 f 3.4 

RGS3 Regulator of G-protein signaling 3 i 2.1 1 2.9 

NUBP1 Nucleotide-binding protein 1 i 2.2 f 2.69 

SLIT3 Slit homolog 3 1 3.3 i21.3 

LOXL3 Lysyl oxidase homolog 3 1 2.73 i 5.7 

CAMPI Secretory carrier- associate membrane protein 1 ^ 3.8 i 3.98 

SEPTS Septin 8 1 2.48 1 3.9 

PPP2R4 Serine/threonine phosphatase 2A, subunit B 1 2.44 1 3.5 

IRS2 Insulin receptor substrate 2 1 2.0 i 3.47 

RAPGEF5 Rap guanine nucleotide exchange factor 5 1 2.01 1 3.3 

ZNFX1 NFXl -type zinc finger containing protein 1 1 2.59 i 3.06 

JPH3 Junctophilin 1 2.0 i 2.6 

TCF4 Transcription factor 4 1 2.04 1 2.31 

IP09 Importin 9 1 2.09 i2.12 

WNK2 Serine/threonine protein kinase WNK2 1 3.6 1 2.1 

Sorting was done by the magnitude of the signal in BA group for both, upregulated ( f ) and down regulated { I ) genes. 
doi:1 0.1 371 /journal.pone.Ol 03740.t004 
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Table 5. Gene ontology annotation analysis: Oxidative stress. 




Symbol 


Gene 


PPA 


BA 


GPX3 


Glutathione Peroxidase, functions in the detoxification of hydrogen peroxide 


t 5.797 


t 7.459 


CYGB 


Cytoglobin, encodes a globin protein found in vertebrate cells 




i 4.574 


APOE 


Apolipoprotein E precursor, essential for the normal catabolism of triglyceride-rich lipoprotein 
constituents 




t 2.28 


IV1T3 


Metallothionein-3, growth inhibitory factor 


t 11.302 


t 8.069 


BNIP3 


BCL2/adenovirus ElB 19 kDa protein-interacting protein 3, protects from virally-induced cell 
death and apoptosis 




I 2.195 


EPHX2 


Epoxide hydrolase 2, binds to specific epoxides and converts them to the corresponding 
dihydrodiols 


t 2.807 


t 4.989 


CYGB 


Cytoglobin 




i 4.574 


DUSP1 


Dual specificity protein phosphatase 1, induced by oxidatlve/heat stress and growth factors, 
may play a role in negativeregulation of cell proliferation 


J, 2.068 


I 2.736 


FOXM1 


Forkhead box protein Ml, encodes transcriptional activator involved in cell proliferation 


J, 5.286 


I 7.152 


NQ01 


NAD{P}H dehydrogenase [quinone] 1, encodes cytoplasmic 2-electron reductase 


J, 3.195 


I 2.749 


Total gene IDs in the ontology- 106, matches for PPA - 9, for BA - 12. 
doi:l 0.1 371/journal.pone.Ol 03740.1005 






Table 6. Gene ontology annotation analysis: Dopamine & serotonin pathways. 








Symbol 


Gene 


PPA 


BA 


COMT 


Catechol 0-methyltransferase catalyzes the transfer of a methyl group from 
S-adenosylmethionine to catecholamines 


I 2.693 


I 2.742 


TH 


Tyrosine 3-monooxygenase, rate limiting enzyme In the synthesis of catecholamines, 
converts tyrosine to dopamine 


i <2 


i 7.069 


DBH 


Dopamine beta-hydroxylase, oxidoreductase which Converts dopamine to norepinephrine 


I 2.157 


i 1 7.434 


DDC 


Aromatic-L-amlno-acid decarboxylase catalyzes the decarboxylation of L-3,4- 
dihydroxyphenylalanine to dopamine, L-5-hydroxytryptophan to serotonin, L-tryptophan 
to tryptamine. 


I 3.181 




MAOA 


Amine oxidase [flavin-containing] A catalyzes the oxidative deamlnation of dopamine, 
norepinephrine, and serotonin 


t 2.049 


t 2.323 


GTPCH 


GTP cyclohydrolase 1, first and rate limiting enzyme in tetrahydroblopterin biosynthesis, 
an essential cofactor for aromatic amino acid hydroxylases and nitric oxide synthases 


t 2.671 




TPH1 


Tryptophan 5-hydroxylase 1 catalyzes the first rate limiting step In the biosynthesis 
of serotonin 


t 2.869 




SLCeA2 


Sodium-dependent noradrenaline transporter 




I 8.594 


ADCY2 


Adenylate cyclase type 2, catalyzes the formation of second messenger cAMP 




t 2.255 


ADCY3 


Adenylate cyclase type 3 




i 2.651 


CASP3 


Caspase-3 precursor, plays a role in the execution-phase of cell apoptosis 




t 2.942 


DUSP1 


Dual specificity protein phosphatase 1 


I 2.068 


I 2.736 


FOS 


Proto-oncogene protein c-fos encodes leucine zipper proteins that dimerize with JUN 
family proteins forming the AP-1 transcription factor complex 


t 5.866 


t 2.22 


PPP1R1B 


Protein phosphatase 1 regulatory subunit IB, encodes a bifunctional signal 
transduction molecule, kinase or phosphatase inhibitor 


t 2.158 


t 54.368 


CYP2D6 


Cytochrome P450 2D6 encodes a member of the cytochrome P450 super family 
of monooxygenases 


t 2.283 


t 6.308 


APP 


Amyloid beta A4 protein precursor encodes a cell surface receptor and transmembrane 
precursor protein cleaved by secretases to form a number of peptides 


i 3.049 


i 5.745 


ARRB2 


Beta-arrestin-2, thought to participate in agonist-mediated desensitization of 
G-protein-coupled receptors and to cause specific dampening of cellular responses to stimuli 


t 2.071 


SLC18A1 


Vesicular monoamine transporter (dopamine gene target} 


i 6.834 


I 82.613 


Total gene IDs in the ontology - 85, matches for PPA - 13, for BA -15. 
doi:l 0.1 371/journal.pone.Ol 03740.1006 
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Table 7. Gene ontology annotation analysis: Fatty acid metabolism. 





Symbol 


Gene 


PPA 


BA 


ACAA2 


3-ketoacyl-CoA thiolase, mitochondrial catalyzes the last step of the mitochondrial 
fatty acid beta-oxidation 


t 2.608 


t 2.982 


AC ATI 


Acetyl-CoA acetyltransferase, mitochondrial, catalyzes the reversible formation 
of acetoacetyl-CoA from two molecules of acetyl-CoA 


t 2.003 




ACAT2 


Acetyl-CoA acetyltransferase, cytosolic encodes enzyme involved in lipid 
metabolism, acetoacetyl-CoA thiolase 


i 4.148 


i 3.199 


ACADSB 


Short/branched chain specific acyl-CoA dehydrogenase, mitochondrial precursor, 
catalyzes the dehydrogenation of acyl-CoA derivatives in the metabolism of 
fatty acids or branch chained amino acids 




t 2.776 


ACADVL 


Very long-chain specific acyl-CoA dehydrogenase, mitochondrial precursor, 
catalyzes the first step of the mitochondrial fatty acid beta-oxidation pathway 




t 2.251 


ACSL5 


Long-chain-fatty-acid-CoA ligase 5, converts free long-chain fatty acids into 
fatty acyl-CoA esters 


i 3.046 


i4.13 


ACSL6 


Long-chain-fatty-acid-CoA ligase 6, catalyzes the formation of acyl-CoA from 
fatty acids, ATP, and CoA 




I 2.283 


ACOT7 


Cytosolic acyl coenzyme A thioester hydrolase, hydrolyzes the CoA thioester of 
palmitoyl-CoA and other long-chain fatty acids 




I 2.274 


CRAT 


Carnitine 0-acetyltransferase, catalyzes the reversible transfer of acyl groups from 
an acyl-CoA thioester to carnitine and regulates the ratio of acylCoA/CoA in the 
subcellular compartments 




I 2.777 


DECR 


1 2,4-dienoyl-CoA reductase, mitochondrial, accessory enzyme in beta-oxidation 
and metabolism of unsaturated fatty enoyl-CoA esters 


t3.7 


■f 3.537 


WlUT 


Methylmalonyl-CoA mutase, mitochondrial, vitamin B12-dependent enzyme which 
catalyzes the isomerization of methylmalonyl-CoA to succinyl-CoA 




i 2.692 


FABP3 


Fatty acid-binding protein, heart, thought to participate in the uptake, intracellular 
metabolism and/or transport of long-chain fatty acids 




i 3.447 


SLC27A1 


Long-chain fatty acid transport protein 1 




t 2.537 


PRKAA1 


5'-AMP-activated protein kinase alpha-1, the catalytic subunit of the 
5'-prime-AMP-activated protein kinase (AlVlPK) a cellular energy sensor 




t 2.715 


BDH2 


3-hydroxybutyrate dehydrogenase type 2 


t 4.612 


t 53.067 


GPD1 


Glycerol-3-phosphate dehydrogenase [NAD+], plays a critical role in carbohydrate 
and lipid metabolism 


J, 2.442 


I 9.572 


LIPE 


Hormone-sensitive lipase 


T 2.238 


t 2.575 


LPL 


Lipoprotein lipase precursor, has dual functions of triglyceride hydrolase and 
ligand/bridging factor for receptor-mediated lipoprotein uptake 




I 8.57 



Total gene IDs in the ontology - 90, matches for PPA - 8, for BA - 18. 
doi:l 0.1 371 /journal.pone.Ol 03740.t007 



including PKA, Ca2-H/CAM kinases and Ras/MAPK pathway 
[78,67]. 

Our data correlate well with the fact that SCFA may enter the 
cell by passive difiusion and/or active transport via specific 
monocarboxylic acid/ketone transporters [79,80] and inhibit 
HDAC, but they can also interact with membrane receptors and 
trigger activation of chfiferent signaling cascades to alter the 
expression of specific genes [81]. In this regard, PPA was 
characterized as the most potent agonist for GPCR41 and 
GPCR43 which are coupled to inositol 1,4,5-trisphosphate 
formation, intracellular Ca2-I- release, ERKl/2 activation 
[82,83]. G-protein signaling coupled to cAMP and regulation of 
adenylate cyclase was also among the top 5 networks identified for 
common genes in our microarray data analyses (see Table 9). In 
rodents and humans, PPA is found to directly regulate GPCR41- 
mediated sympathetic outflow and thereby to modulate body 
energy expenditure in maintaining metabolic homeostasis [84,85] . 



PPA nnodulation of inflammatory processes/lipid 
homeostasis 

PPA can also exert its effect by activating nuclear receptor 
PPAR gamma [86]. We have demonstrated previously that PPAR 
gamma receptors are expressed in PC 1 2 cells and are involved in 
the independent regulation of catecholaminergic and opioid 
pathways by SCFA and related drugs like 4-phenylbutyrate [75]. 
Our microarray analyses showed SCFA - altered expression of 
genes, involved in fatty acid metabolism and transport, ketogenesis 
and ketone body metabolism, as well as triacylglycerol metabolism 
(Table 7). These findings are of particular interest since abnormal 
lipid metabolism has been reported in ASDs and ASD-related 
disorders [87,88,89,90,91,92]. Moreover, clinical studies suggest 
improvements in core symptoms of patients with ASDs following 
supplementation with polyunsaturated fatty acids (PUFAs) 
[93,94,95] and we observed altered lipid profiles in rat brain 
phospholipids following infusion with both PPA and, to a lesser 
extent, BA [37,96]. 
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Table 8. Gene ontology annotation analysis: GAP junction proteins. 





Symbol 


Gene 


PPA 


BA 


GJB5 


Gap junction beta-5 protein, involved in intercellular communication related to 
epidermal differentiation and environmental sensing 




i 2.624 


EGFR 


Epidermal growth factor receptor precursor 




i 2.688 


PDGFRA 


Alpha-type platelet-derived growth factor receptor 




i 2.152 


DBN1 


Drebrin, cytoplasmic actin-binding protein thought to play a role in the process 
of neuronal growth 


f 2.128 




CTNNB1 


Catenin beta-1, part of a complex of proteins that constitute adherens junctions 




t 2.21 


TJP1 


Tight junction protein ZO-1 




t 2.687 


PLCB1 


1-phosphatidylinositol-4,5-bisphosphate phosphodiesterase, catalyzes the formation 

of Inositol 1,4,5-trisphosphate and diacylglycerol from phosphatidylinositol 4,5-bisphosphate 




t 2.663 


SRC 


Proto-oncogene tyrosine-protein kinase Src, involved in regulation of embryonic 
development and cell growth 


i 2.286 


i 5.815 


TUBB2B 


Tubulin beta-2B chain, binds GTP and is a major component of microtubules. 


1 3.164 




TUBB2C 


Tubulin beta-2C chain 


i 2.047 




GRB2 


Growth factor receptor-bound protein 2, binds the epidermal growth factor receptor, 
involved in signal transduction 




i 3.034 


DAC1 
KMr 1 


RAF proto-oncogene serine/threonine-protein kinase, encoded protein is a MAP 
kinase kinase kinase (MAP3K} 




J, 2.462 


ITPR2 


Inositol 1,4,5-trisphosphate receptor type 2 


^4.719 


t 8.015 


PLCB1 


1-phosphatidylinositol-4,5-bisphosphate phosphodiesterase 




t 2.663 


PRKCG 


Protein kinase C gamma type, expressed solely in the brain and spinal cord, 
localization is restricted to neurons 




T 2.203 


ADCY2 


Adenylate cyclase type 2 




t 2.255 


ADCY3 


Adenylate cyclase type 3 




i 2.651 



Total gene IDs in the ontology - 101, matches for PPA - 6, for BA- 12. 
doi:1 0.1 371 /journal.pone.Ol 03740.t008 



PPA activation of CREB dependent catecholaminergic 
pathways 

The microarray data reported here are also consistent with the 
activation of CREB pathway - i.e. increased expression of 
CREB3LI was observed in PPA (2.33 fold) and BA (4.2 fold) 
group, accompanied by elevated MAP4K1 (2-fold and 3.6 fold 
resp.), CAMK2G (2-fold for both SCFA) and CAMKV (increased 
by approximately 3-fold by PPA and BA). Increased P-CREB 
immunoreactivity was observed in ASD- hippocampus and white 
matter regions following intracerebroventricular infusion of SCFA 
(ASD rodent model, [35,3]). Moreover, our prior work with BA is 
consistent with a signal transduction system involving CREB and 



converging on the CRE and a novel upstream 5' regulatory 
element (BRE; GCCTGG at -509 to -504 of both, the rat TH 
and proenkephalin promoters [60,58,55,56] . 

The ability of PPA to activate endogenous TH transcription at 
low concentrations ( 1 mM) results in accumulation of TH mRNA 
and TH protein (Fig. 2) indicative of increased cell capacity to 
produce catecholamines. This observation is similar to the effects 
of BA and VPA in PC 12 cells we reported before [63,58]. 
Furthermore, it provides a possible molecular mechanism to 
account for the hyperactivation of the mesocortical dopamine 
system suggested in several clinical reports of ASD [97,98] and in 
animal models of ASD [66,99]. Other 4 carbon structures (i.e. 
GABA or y-OH butyrate - which is accumulated as a result of 



Table 9. Most relevant biological networks identified for common genes. 



# Processes Pathways P-value Zscore 



1 


Cell cycle, mitotic cell cycle 


230 


1.73e-19 


15.25 


2 


Cell cycle, regulation of progression through cell cycle 


147 


1.16e-15 


13.07 


3 


Intracellular signaling cascade, localization of the cell, cell motility 


71 


1 .84e-5 


5.90 


4 


G-protein signaling coupled to cAMP, regulation of adenylate cyclase 


52 


4.02e-17 


13.97 


5 


Actin filament polarization and/or depolarization, protein polarization 


0 


2.98e-68 


34.11 



The list of biological networks was generated using the Analyze Networks (AN) algorithm with default settings (MetaCore), based on relative enrichment with the 
uploaded data and the relative saturation of networks with canonical pathways. 

P-value: represents the probability for a particular mapping of an experiment to a map {or network, or process) to arise by chance considering the number of genes in 
the experiment versus the number of genes in the map within the "full set" of all genes on maps. Z-score ranks the sub-networks of the AN algorithm with regard to 
their saturation with genes from the experiment. 
doi:1 0.1 371/journal.pone.01 03740.t009 
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1 . Synaptic transmission 

2. Regulation of GPCR protein signaling 

3. Response to stimulus 

4. Transmission of nerve impulse 

5. Cardiac muscle development 

6. Regulation of norepinephrine secretion 

7. Norepinephrine secretion 

8. Negative regulation of caspase activity 
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Common genes 
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Figure 4. Enrichment analysis of differentially expressed genes: distribution by gene ontology (GO) processes. Differentially expressed 
genes in BA and PPA groups (t-test compared to control group p<0.01 ) were subjected to enrichment analysis (which consists of matching gene IDs 
for the common, similar and unique sets of the uploaded files with gene IDs in functional ontologies in MetaCore). The figure illustrates the 
distribution by GO processes. The gene content is aligned between all uploaded files. The set of common gene IDs is marked as blue and white 
stripes. The unique genes for the files are marked as colored bars (BA - orange; PPA- blue). The sorting was done by common gene IDs; p-value was 
set for 0.05; both signals (induced and repressed) were included. The data shown are for sorting method "similarity by". The degree of "relevance" to 
different categories for the uploaded datasets is defined by p values, so that the lower p-value gets higher priority. The top 10 processes are listed 
based on their -log (p-value). 
doi:10.1371/journal.pone.0103740.g004 



succinic semialdehyde dehydrogenase deficiency, a rare autosomal 
recessive disorder which presents with autistic beliaviors, [100]) 
and following ketogenic diet administration, a possible therapy for 
ASD, also have potent effects on CNS function via the 
dopaminergic system suggesting a broader implication of this 
putative mechanism. 

Although the stimulatory effect of SCFA on the TH promoter is 
proportional to the drug dose used (Fig. lA for PPA, also 
previously shown by us for BA, phenylbutyrate and VPA 
[75,58]), TH mRNA and protein fall below basal levels at high 
concentration [58,63], which was associated with increased TH 
mRNA destabUization presumably via a butyrate response factor 
(BRF) acting at a 3' untranslated region [58,61,63,101]. These 



high concentrations are comparable to the doses of brief pulse 
intracerebroventricular infusions used in our rat ASD model [35] 
and blood PPA concentrations in human propionic acidemia 
(4400 uM/L, see [25,27]). An interesting new finding of the 
present study is the bi-directional perturbation of the dopami- 
nergic pathways by both PPA and BA (Table 1, 6 and Table S3). 
It should be emphasized that this novel SCFA-triggered mode of 
regulation is different from the reported multiple studies (rev. in 
[102,103,104,59]) where diversified stimuli (cyclic AMP analogues, 
phorbol esters, growth hormones, glucocorticoids, depolarization 
agents) each increased TH mRNA levels and stimulated TH 
promoter activity in PC 12 cells. Whether a similar bi-directional 
dysregialation of the brain catecholaminergic system occurs in 
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Table 10. SCFA alter the expression of genes implicated in ASD in human and animal studies. 





Genes implicated in ASD (Copy number variations) 


Symbol 


Gene 


PPA 


BA 


C3orf58 


Uncharacterized protein C3orf58, renamed deleted in autism 1, encodes signal 
peptide for targeting to the secretory pathway 




t 4.539 


CACNA1C 


Voltage-dependent L-type calcium channel alpha-IC mediate the influx of calcium 
ions into the cell upon membrane polarization 




i 2.093 


DMPK 


Disks large-associated protein 2, dystrophy myotonica protein kinase 


f 2.258 


] 3.018 


FMR1 


Fragile X mental retardation 1 protein, negative regulator of translation, plays 
key role in neuroplasticity 


i 2.289 


i 2.383 


NIPBL 


Nipped-B-like protein, facilitates long distance enhancer-promoter communications 


i 2.766 




NLGN2 


Neuroligin-2 , post synaptic cell adhesion molecule of inhibitory synapses 




i 2.471 


NRXN1 


Neurexin-1 -alpha, 1-beta, pre-synaptic protein that forms complexes with 
neuroligins or other partners to shape synaptic neurotransmission 


i 2.227 




PTEN 


Phosphatidylinositol-3,4,5-trisphosphate 3-hosphatase and dual- specificity 
protein phosphatase, tumor suppressor negatively regulate AKT/PKB signaling 




t 2.922 


RAI1 


Retinoic acid-induced protein 1, transcription factor induced by retinoic acid 




i 2.566 


ASD Candidate genes: 


ASTN1 


Astrotactin-1 precursor, neuronal adhesion molecule, involved in glia-guided neuronal 
migration 




i 2.68 


DLGAP2 


Disks large-associated protein 2, neuronal postsynaptic density protein 


t 2.258 


t 3.018 


DPP6 


Dipeptidyl aminopeptidase-like protein 6 auxiliary subunit of Kv4-mediated A-type K(-i-} 
channels 


i 2.455 


i 1 5.024 


DPYD 


Dihydropyrimldine dehydrogenase [NADP-h], involved in uracil and thymidine catabolism 




t 5.115 


GABRB1 


Gamma-aminobutyric acid receptor subunit beta-1 




t 2.368 


KCNMA1 


Calcium-activated potassium channel alpha-1 the channels control neuronal 
excitability and smooth muscle tone 


i 2.885 




MAOA 


Monoamine oxidase [flavin-containing] A degrades amine neurotransmitters 


t 2.049 


f 2.323 


MTFI 


Metal regulatory transcription factor 1, Zn binding transcription factor that mediates 
regulation of metallothionein genes 




t 2.056 


PCDH9 


Protocadherin-9, encodes a cadherin-related neuronal receptor that localizes to synaptic 
junctions 




] 22.029 


PCDH10 


Protocadherin alpha 10, cadherin-related neuronal receptor believed to play a role in 
establishment and function of specific cell-cell connections in the brain 




i 2.477 


PIP5K1B 


Phosphatidylinositol-4-phosphate 5-kinase 1 beta 


t 2.328 




PRKAG2 


5'-AMP-activated protein kinase gamma-2, energy-sensing enzyme, composed 
of a catalytic alpha, noncatalytic beta and noncatalytic regulatory gamma subunit. 


i 13.619 


i 2.1 54 



doi:l 0.1 371/journal.pone.Ol 03740.t01 0 

response to excessive concentrations of SCFA exist in vivo remains 
to be examined. A bi-directonal system is intriguing since high 
concentrations of SCFA down regulated not only TH, but also 
DBH, DDC and COMT genes in PC 12 cells (our results reported 
here, also [62]). Loss or down regulation of the DBH gene is 
strongly associated with deficits in affiliative and self injurious 
behaviors in several animal models (rev. in [105,106]) and with 
ASD phenotype in humans [107,108,109,110]. 

Serotonergic/Cholinergic systems 

Of interest, only PPA was able to alter the serotonin system in 
PC 1 2 cells by inducing the expression of Tryptophan 5-hydoxylase 
1 (TPHl, rate-limiting enzyme in the synthesis of the neurotrans- 
mitter serotonin) and GTP cyclohydrolase 1 (GTPCH - the rate 
limiting enzyme in BH4 biosynthesis, an essential cofactor of three 
aromatic amino acid hydrolases including TPH and TH); Table 6. 
Dysfunctional serotonin signaling has been suggested as an 
etiologic factor in ASD (rev. in [111]). Serotonin is linked to the 
mediation of several psychological processes, many of which are 



altered in ASD patients, including mood, social interaction, sleep, 
obsessive compulsive behaviors and aggression. Similar symptoms 
were also observed in our PPA animal model of ASD [3,44]. 

During the past few years, studies have begun to suggest that 
cholinergic systems in the brain might also be implicated in ASD 
[112,113]. We also found altered expression of CHRNA3, 
CHRNBl (both down regulated by BA), CHRNA5 and CHRNB4 
(up-regulated by PPA 3.2 and 5-fold, see Tables SI and S2, and 
GSE56516). 

Mitochondrial dysfunction/oxidative stress 

The brain is vulnerable to mitochondrial dysfunction, oxidative 
stress and excitotoxicity [1 14]. Consistent with the effects of SCFA 
(PPA and/or BA) on ceUular metabolism [35,37,96,1 15], the data 
presented in Table 5 are similar to the clinical presentation of 
patients with ASD who show metabolic dysfunction, including 
impairments in B12, glutathione, or carnitine metabolism 
[116,117,118,43] and mitochondrial disorder/dysfunction 
[119,120]. VPA similarly alters mitochondrial metabolism and 
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causes depletion of carnitine stores and encephalopathy [121,122]. 
Of note, we have shown specific alterations of carnitine and 
mitochondrial metabolism in a large cohort of ASD patients 
predicted with the PPA rodent model suggesting mitochondrial 
dysfunction, which also alters gene expression [123], may play a 
key role in ASD pathogenesis [45,3,4]. 

General effects of SCFA 

DNA microarray technolog}' enables simultaneous detection of 
perturbations in thousands of genes in a single experiment. The 
present study also identified a large number of affected genes 
following exposure of PC 12 cells to high concentrations of SCFA 
using a criteria of S 2-fold change and a PsO.05 (Figure 3, Tables 
SI, S2, S3): 1599 unique for BA, 355 unique for PPA and 1010 
common for both SCFA. Upon closer look at the expression 
profiles triggered by PPA and BA, we found a set of 16 common 
genes altered to a higher extent after PPA treatment (Table 3) and 
a total of 2 1 common genes affected in the opposite direction by 
PPA and BA (Table 4). Together, these data confirm the notion 
that in addition to the common ability to cause HDAC inhibition, 
other proct'sscs in PC 1 2 cells are influenced in a PPA-and/ or BA 
specific marmer providing further layers of regulation and 
potential for neuronal plasticity. 

Alteration of ASD associated genes FMR1, Neurexins, 
Neuroligins 

In addition, we identified a number of genes (listed in Table 1 0) 
that were previously implicated in ASD or candidate genes for 
ASD from evaluation of copy number variation genetic studies. 
For example, loss of function of Fragile X Mental Retardation 1 
gene (FMRl) causes Fragile X syndrome with up to 90% of 
affected children exhibiting ASD symptoms [124,125,126]. 
Treatment of PC 1 2 cells with PPA or BA caused down regulation 
of FMRl gene expression (Table 10). Neurexin 1 (NRXNl) is 
another gene considered to be causal for ASD [127,128] which 
was also down regulated following BA administration. Many genes 
associated with ASD are also involved in the neuroligin-neurexin 
interaction at the glutamate synapse [129,28]. 

Thus, it appears that number PPA and BA affected genes are 
related to ASD. In other words, the apparent epigenetic regulation 
of gene expression by gut-microbiota products (SCFA) can in 
principle result in similar "loss" or "gain" of function and disrupt 
similar pathways/networks as genetic alterations associated with 
ASD. 

We applied MetaCore Ceiie Ontology' (GO) bioinformatics tool 
[130] for functional enrichment profiling of our genome-wide 
differential expression data in response to SCFA (Figure 4). Of 
note, the top 10 enriched GO terms included synaptic transmis- 
sion, regulation of G-protein coupled signaling, transmission of 
ner\'e impulses. Our data are concordant with the comprehensive 
review and analysis of published literature and data of ASD 
genetics [131] including a total of 2193 genes, 2806 SNPs/ 
VNTRs, 4544 CNVs and 158 linkage regions. 

In summary, the gene expression analysis in this study provides 
valuable insight into the complexity of possibilities and specific 
molecular pathways or processes altered in PC 12 cells upon 
stimulation with SCFA, many of which, directly or indirecdy, have 
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capacity of SCFA to regulate multiple pathways and networks, 
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that our data represent a useful initial cell model resource for 
identifying and functionally testing SCFA-responsive neuronal 
genes. 

Although the preliminary data are compelling, it is important to 
note that the effects of SCFAs on cellular systems are complex, 
often indirect, and often time, tissue- and dose-related [8,3,132]. 
In addition to direct HDACI effects, other SCFA effects on 
GPCR, second messenger/ calcium signaling, and mitochondrial/ 
oxidative stress also alter gene expression. One should be cautious 
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research include using whole animal models at specific time 
periods and routes of PPA exposure to examine the effects of gut 
derived bacterial metabolites on brain function and behavior, and 
their possible relation to neurodevelopmental disorders. 

We speculate that altered levels of environmental factors (such 
as SCFA) normally produced in the gut by microbial fermentation 
of dietary' carbohydrates, may: 

1) Represent a novel evolutionary link to neural plasticity 
associated with different behavioral responses to a changing 

environment. 

2) Affect brain and behavior by multiple molecular mechanisms, 
some of which involve Ca2-I-/ cAMP/ CREB signaling path- 
ways. 

3) Modulate expression of key genes (neurotransmitters, cell 
adhesion molecules, genes involved in lipid metabolism, 
oxidative stress and mitochondrial function) many of which 
have been implicated in neurodevelopmental disorders, such 
as ASD. 
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